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ABSTRACT 

The starbusting, nearby (D = 32.9 Mpc) spiral (Sc) galaxy NGC 2276 belongs to the 
sparse group dominated by the elliptical galaxy NGC 2300. NGC 2276 is a remarkable 
galaxy, as it displays a disturbed morphology at many wavelengths. This is possibly due 
to gravitational interaction with the central elliptical galaxy of the group. Previous 
ROSAT and XMM-Newton observations resulted in the detection of extended hot 
gas emission and of a single very bright (^10 41 erg s -1 ) ultraluminous X-ray source 
(ULX) candidate. Here we report on a study of the X-ray sources of NGC 2276 based 
on Chandra data taken in 2004. Chandra was able to resolve 16 sources, 8 of which 
are ULXs, and to reveal that the previous ULX candidate is actually composed of a 
few distinct objects. We construct the luminosity function of NGC 2276, which can be 
interpreted as dominated by high mass X-ray binaries, and estimate the star formation 
rate (SFR) to be ^5-15 Mq yr -1 , consistent with the values derived from optical 
and infrared observations. By means of numerical simulations, we show that both ram 
pressure and viscous transfer effects are necessary to produce the distorted morphology 
and the high SFR observed in NGC 2276, while tidal interaction have a marginal effect. 

Key words: galaxies: individual: NGC 2276 - galaxies: star formation - X-rays: 
binaries - X-rays: galaxies - methods: numerical. 


1 INTRODUCTION 

NGC 2276 is an Sc galaxy belonging to the poor group 
of galaxies NGC 2300, which is composed of four or five 
galaxieQand dominated by the namesake elliptical galaxy. 
Here and throughout the paper we assume a distanc^] to 
NGC 2276 of 32.9 Mpc and rescale published values to this 
distance when necessary. The NGC 2300 group is the first 
in which an X-ray emitting intragroup medium (IGM) was 
observed ( Mulchaey et al.| 1993). Observations with ROSAT 


* E-mail: anna.wolter@brera.inaf.it 

1 Apart from NGC 2300 and NGC 2276, the other members are 
NGC 2268 and IC455 flHuchra fc Geller|1982| . UCG 03670, at 16' 
(~0.15 Mpc) from NGC 2300, is likely another group member. 

2 From the NASA/IPAC Extragalactic Database 
(NED) http://ned.ipac.caltech.edu/ and adopting Ho = 
73 km s -1 Mpc -1 


in fact revealed an extended (c^0.2 Mpc), unusually dense 
5.3 X 10 -4 cm -3 ), hot (~0.9 keV) and relatively metal 


poor (~0.06 Zq) intragroup gas halo ( Mulchaey et al.| 1993 


Davis et al.|19 96). NGC 2276 is undergoing starburst activ¬ 


ity. A number of Hu regions were identified in the galaxy, 


spread throughout its volume (Hodge &; Kennicutt 1983 
Davis et al. 1997), and many super novae have been observed 


the last fifty years (e.g. Iskudaryan & Shakhbazyan 1967 


Barbon, Cappellaro &; Turatto||19891 |Treffers et al. 1993 
Dimai, Migli ardi &; Man zini||2005|). The star formation rate 


(1983 ) to be 9.5 Mq yr 


-l 


(SFR) was estimated by Kennicutt ( 
based on the high Ha luminosity, and values between 5 and 
15 Mq yr -1 can be found in the literature from different 
measurements (mostly Ha and infrared observations; e.g. 
Ken nicutt Jr. , Tambl yn &; Con gdon 1994, Sand ers et al.| 
2003||James et al.||2004|). 

The asymmetric shape of NGC 2276, with a ‘bow-shock- 
like’ structure along its western edge and a ‘tail’ of gas ex- 
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tending towards east, and the unusually intense star forma¬ 
tion, especially along the western edge, (the source falls in 
the upper range of the distribution of star-forming galaxies; 
see e.g. |Noeske et al.|2007| for the scatter-plots of galaxies in 
the Aegis field) have attracted considerable attention but, 
while it is quite clear that these features are due to inter¬ 
action with the environment, no consensus on the dominant 
physical process involved has yet been reached. 

Since NGC2300 and NGC2276 are relatively close to 
each other (7', about 70 kpc), a prime candidate is the tidal 


The Chandra data allowed Rasmussen et al. (2006) to 


interaction between the two galaxies (Gruendl et al. 1993 
Davis et al.||1996~ 1997 Elmegreen et al.||1991 ). Other pro¬ 

cesses involved in shaping NGC2276 are the viscous and 
the ram-pressure stripping of the galactic gas in the IGM 


(Nulsen|1982 

Rasmussen, Ponman X Mulchaey|2006). Ras- 

mussen et al. 

(2006) questioned the importance of the tidal 


interactions and argued that turbulent viscous and ram- 
pressure stripping, ram pressure compression of the disc gas 
and starburst outflows are more likely the main elements 
influencing the shape and SFR of NGC 2276. 

Using an XMM-Newton observation carried out in 
2001, Davis X Mushotzky ( |2004| identified a bright X- 
ray source (XMMU J072649.2+854555) at the western edge 
of the galaxy with an ultraluminous X-ray source (ULX: 
non-nuclear sources with X- ray lum i nosities in excess of 


L x = 10 ' 


39 erg s 1 


see e.g. 


Feng X Soria 


2011 


for a re¬ 


view) candidate. The X-ray luminosity of this object in 


the 0.5-10 keV band was L 


0.5-10 keV 


= 1.1 X 10 £5 7 = 


5.7 x 10 4 O d | 2 .9 erg s -1 (where d/v stands for the distance 
in units of N Mpc), making it one of the most luminous 
known in its class (e.g. [Sutton et al. 2012 ). They also ob¬ 
served that the nuclear source, undetected in the XMM- 
Newton data, was much dimmer (by a factor of ~ 10 at least) 


than the ROSAT /HRI nuclear source seen by Davis et al. 
(1997) eight years earlier at L^ 5 ~ 2keV ~ 2 x 10 4 U dJ”^ = 
10 4U d32.9 erg s -1 . 

ULXs are mainly found in galaxies with high SFRs, and 
indeed most might be assumed to be the extension of the 
High Mass X-ray Binaries (HMBX) population to high X-ray 


luminosities. For instance, in The Antennae (Zezas X Fab- 
biano ||2002 ), in the Cartwheel galaxy (Wolter X Trinchieri 


2004) , in the NG C 2207/IC 2163 pair ( [Mineo et ahp 013 ) 
and in NGC 4088 ( Mezcua et al.|[2014 ) the X-ray luminos¬ 


ity function follows the “universal” law found by |Grimm| 


et al. (2003) for HMXBs when scaled for the SFR. This 


assumption, linked to the proposed influence of low metal- 
licity on the evolution of massive stars (see |Mapelli et al. 


2010 and references therein) leads to expect not one, but 


a large number of ULXs in NGC 2276. Using the corre¬ 
lations by Mapelli et al.| (2010J and assuming a metallic- 


ity Z m 0.22 Zq and a SFR of 10 Mq yr _ , we expect 
about 10 ULXs. The metallicity is derived applying the 
|Pilyugin X Thuan] ( |2005) calibration on the de-reddened 
[Oil A 3727]/H/3 and [OIII A5007]/H/3 ratios measured by 
Kennicutt ( |1992) on an integrated spectrum of NGC2276. 
Indeed, a high-spatial-resolution observation of NGC 2276 
taken with Chandra (Obs.ID 4968, PI: Rasmussen) revealed 
that the bright ULX candidate is actually composed of a few 
different sources ( Wolter et al.|201l] ), and showed a total of 
8-9 sources bright enough (Lx > 10 39 erg s -1 ) during the 
observation to be classified as ULXs (Wolter et al.|2011 Liu 

WI). 


resolve details of the X-ray morphology of NGC 2276, which, 
as seen at other wavelengths, is characterised by a shock-like 
feature along the western edge and a low surface brightness 
tail extending to the east. Spatially resolved spectroscopy 
shows that the data are consistent with the disc gas being 
pressurized at the leading western edge of NGC 2276, due 
to the galaxy moving super sonically through the IG M, at 
a velocity of ~900 km s - 


Rasmussen et al. 


(2006) esti¬ 


mated that the diffuse hot gas in the disc of NGC 2276 has 
a temperature of kT ~ 0.3 keV and an X-ray luminosity 
of £^- 3-2 - 0keV ) — 1.9 x 10 4 O d § 6 . 8 = 1.5 x 10 4 O di 2 . 9 erg s -1 , 
with an inferred residual component of unresolved X-ray bi¬ 
naries of about 10-15 per cent of the total flux. The IGM, 
as observed by Chandra , has a profile consistent with the 
previous ROSAT measurements ( Davis et al.||1996 ), a tem¬ 
perature of kT « 0.9 keV and Z = 0.17 Zq (Rasmussen 


et al.|2006 ). This metallicity is well matched, given also the 
different method used, to the one used above 


In Wolter et al. (2011), we reported on the several new 


ULXs observed by Chandra in NGC 2276 and on the fact 
that the bright ULX candidate XMMU J072649.2+854555 
actually consists of a number of distinct sources (see also 
Sut ton et al.||2012 ). In this work, we give more details on the 
point sources, presenting the whole Chandra data analysis, 
investigate the main mechanisms at work in shaping the 
galaxy by means of numerical simulations, and discuss their 
influence on the X-ray source population. 


2 ANALYSIS OF CHANDRA DATA 


The Chandra observation (Obs. ID 4968, see Rasmussen] 
et al.||2006 for further details) was performed on 2004 June 
23 with the Advanced CCD Imaging Spectrometer (ACIS; 
Garmire et al. 2003), and lasted about 45.8 ks. The ACIS 
was operated in the standard timed exposure full-frame 
mode, with the ‘very faint’ telemetry format, and NGC 2276 
was positioned on the back-illuminated chip S3 (CCD7). The 
data were reprocessed with the Chandra Interactive Analysis 
of Observations software package (CIAO, version 4.6). 0 We 
removed the pixel randomization, since it deteriorates the 
spatial resolution, but, other than that, we followed stan¬ 
dard data reduction and analysis procedures. 

The first step of our analysis was to run the CIAO 
wavelet source-detection algorithm with SlGTHRESl]^] set to 
10 -6 , so that we expect at most one spurious source in the 
area under examination given that we consider only a single 
CCD. The routine found 41 sources in CCD7, of which 16 


in the galaxy area, defined by D 25 = 177. 1 (de Vaucouleurs 
et al.||1991 ) with a signal-to-noise ratio higher than 4. 


Source positions (in RA order) and statistics of the 16 
sources are given in Table[l] while in Fig. [l] they are shown in 
the X-ray image and over-plotted on an optical image of the 
galaxy. Source counts were extracted from circular regions 
centered at the positions found by the source-detection rou¬ 
tine and with radius 2 ". The sources are all at < 2' from the 
aimpoint, therefore the 2" circle is large enough to include 


3 http://cxc.harvard.edu/ciao/ 

4 See http://cxc.harvard.edu/ciao/ahelp/wavdetect.html for the 
parameters of wavdetect. 
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Figure 1 . Left: Smoothed image (csmooth) from Chandra data in the 0.5-7 keV band. The green circle (radius=50") indicates the 
region used to extract the diffuse component spectrum (associated to S7, see text for details) Right: The X-ray sources’ positions are 
over-plotted on the optical image of NGC 2276 taken from the Digitized Sky Survey (DSS), see http://archive.stsci.edu/dss/index.html 
Open circles and numbers identify detected point sources; red sources are ULXs. White sources have L^° _1U ^ 10 39 erg s -1 . 


Table 1 . Source detected in NGC 2276. Positions from the de¬ 
tection algorithm and net counts from the extraction regions of 
2" radius (see text). No astrometric correction has been applied 
and therefore the positional uncertainty is the standard Chandra 
one of A ^ 0'/6 (90% c.l.). 


Number RA (J2000) 

Dec (J2000) 

Net Counts 
(0.3-5 keV) 

SI 

07 h 26 m 37?4 

+85°45 , 34'/5 

417.4 ±20.4 

S2 

07 h 26 m 43?5 

+85°45'007l 

46.2 ±6.9 

S3 

07 h 26 m 46?l 

+85 0 45'38'/5 

43.2 ±6.6 

S4 

07 h 26 m 47?9 

+85°45'52'. , 3 

175.8 ± 13.3 

S5 

07 h 26 m 48?l 

+85°45'54 , /9 

338.6 ± 18.4 

S6 

07 h 26 m 48?2 

+85°45'49'.'0 

104.7 ± 10.2 

S7 a 

07 h 27 m 13?0 

+85°45'16 , /0 

40.3 ±6.4 

S8 

07 h 27 m 14?9 

+85°46'11'. , 5 

14.2 ±3.9 

S9 

07 h 27 m 15?4 

+85°45'54 , /9 

8.2 ±3.0 

S10 

07 h 27 m 19?7 

+85°46'32'. , 9 

19.2 ±4.5 

Sll 

07 h 27 m 25?8 

+85°45'25 , /9 

68.2 ±8.3 

S12 

07 h 27 m 28?7 

+85°45 , 11'. , 1 

21.3 ±4.7 

S13 

07 h 27 m 53?3 

+85°46'08'/8 

53.2 ±7.3 

S14 

07 h 27 m 58?4 

+85°44 , 37'. , 9 

19.3 ±4.5 

S15 

07 h 28 m 16?0 

+85°44'36'. , 4 

29.2 ±5.5 

S16 

07 h 28 m 19?7 

+85°44 , 28'. , 5 

13.3 ±3.7 


a Source S7 is the center of the galaxy: for subsequent analysis 
counts have been extracted in a circle of 50" radius, excluding all 
detected sources, resulting in a total of 1413 d= 67 net counts. 


more than 90% of all counts [^] and is a very good match to 


the ellipse from the detection algorithm. We expect a very 
small contribution in the area from both the instrumental 
and the diffuse background. The dominant source of error 
is statistic. For source S7, which corresponds to the nuclear 
region, we extracted counts in a circle of 50" radius (ex¬ 
cluding all the detected sources in that area) for consistency 
with the Rasmussen et al. (2006) analysis (see Section 4). 
The sources’ positions are plotted in Fig.[l]and the extrac¬ 
tion area for S7 is indicated as the green circle. Instrumental 
and cosmic background was estimated from circular regions 
outside of the galaxy, but as close as possible to it, and ex¬ 
cluding detected sources. We did not attempt to subtract 
the diffuse gas contribution, which is, at any rate, negligi¬ 
ble in the extraction area of the point sources (less than 1 
count per source). A first indication of the sources’ spectral 
shape can be evinced from Fig. [2] in which the “true” color 
image of the X-ray emission has been constructed. The ab¬ 
sorbed or flat sources appear as blue/green in this image.For 
a more formal analysis, we extracted individual spectra for 
the sources with more than a hundred net counts (SI, S4, 
S5, and S6, plus the extended region around S7). For the 
other sources (S2, S3, S8, S9, S10, Sll, S12, S13, S14, S15, 
S16), since the paucity of counts hinders individual fits, we 
extracted a cumulative spectrum, which contains 336 ± 20 
net counts. All the spectral analysis was performed with the 
xspec (v. 12.8.0) fitting package (Arnaud 1996). Results and 
analysis of the spectra are described in the next Sections. 


5 See e.g. fig. 4.7 in http://cxc.harvard.edu/proposer/POG/html/chap4.html 
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Table 2. Results of spectral fit for the 4 ULXs with more than 100 net counts and for the sum of the fainter 
sources. Unabsorbed fluxes and luminosities are in the 0.5-10 keV energy band. 


Name 

r 

N H 

(10 21 cm” 2 ) 

£-,(0.5 — 10 keV) 
r x 

(IQ- 14 er g cm -2 s -l) 

£(0.5 — 10 keV) 
(10^ 9 erg s -1 ) 

£(2-10 keV) 
(10 3 Urg s -1 ) 

xl (dof) 

SI 

1.96 ±0.34 

4.1 ±2.0 

14.45 ±0.71 

18.7 

10.29 

0.92 (17) 

S4 

1 48+ 0 - 49 
± -^°-0.25 

< 25.0 

6.78 ±0.51 

8.73 

6.39 

1.47 (5) 

S5 

1.77 ±0.40 

3.1 ±2.0 

11.10 ±0.60 

14.3 

8.90 

1.54 (13) 

S6 a 

> 3.2 

6.0 ±4.0 

2.60 ±0.25 

3.36 

2.86 

1.50 (1) 

SUM 

1.93 ±0.35 

1.7 ±0.9 

6.00 ±0.37 

7.74 

4.41 

1.11 (18) 


a Flux and luminosities were computed with T = 3.2 and the fitted value for Nn= 25.3 X 10 21 cm 2 . 



Figure 2. Three color image of the Chandra data. Smoothed with 
CSMOOTH from ciao in the 0.3-0.9 keV, 0.9-1.5 keV and 1.5-8 keV 
bands for red, blue and green respectively. Many point sources ap¬ 
pear blu/green indicating excess absorption and/or hard spectra. 


3 POINT SOURCES AND ULXS 


The distribution of the detected sources follows roughly the 
position of the arms of NGC 2276 and the zones of high ac¬ 
tivity in the other bands. Contamination from background 
sources may be estimated according to the log JV—log S dis¬ 


tribution of Hasinger et al. (1993). At the detection limit of 

_ —2 -1 


erg cm “ s ^ (which corresponds to 


^0.6-2 keV) = 7xl0 —16 

Lx « 10 38 erg s -1 at the distance of NGC 2276) we expect 
1.4 sources by chance in the total area covered by the galaxy. 
Therefore, at least one of the fainter sources is probably a 
background source. 

The few bright sources with enough counts (~ 100-400) 
to allow individual spectral fits are modeled by a simple 
power law modified for the interstellar absorption. A satis¬ 
factory description of their spectra has a hydrogen column 
density Nu of a few 10 21 cm -2 (the total Gal actic Ng in t he 
direction of NGC 2276 is ~5.7 x 10 20 cm -2 ; Kalberla et al. 


2005) and slopes between T = 1.4 and 2.2. The observed nor¬ 


malised net count distribution with the best fit model is plot¬ 
ted in Fig. [3] while the fit results are listed in Table[2] These 
spectra are consistent with typical results for ULXs when 
only low-count-statistics spectra are available (e.g. Swartz 


Table 3. Point source unabsorbed fluxes and luminosities com¬ 
puted from the power-law fit to the mean spectrum of the fainter 
sources (see Sections [ 2 ] and [3j . 


Name 

£-,(0.5 — 10 keV) 
(IQ- 14 er g cm -2 s -l) 

T (0.5-10 keV) 

(IQ 39 er g g-1) 

T (2-10 keV) 

(10 39 erg s -1 ) 

S2 

0.96 ±0.14 

1.24 

0.70 

S3 

0.90 ±0.14 

1.16 

0.66 

S8 

0.30 ±0.08 

0.38 

0.22 

S9 

0.17 ±0.06 

0.22 

0.13 

S10 

0.40 ±0.09 

0.52 

0.29 

Sll 

1.42 ±0.17 

1.83 

1.04 

S12 

0.44 ±0.10 

0.57 

0.32 

S13 

1.11 ±0.15 

1.43 

0.81 

S14 

0.40 ±0.09 

0.52 

0.29 

S15 

0.61 ±0.11 

0.79 

0.44 

S16 

0.28 ±0.08 

0.36 

0.20 


et al. 2004). For the other sources, we considered the cumula¬ 
tive spectrum (see Section|2j for a total of ~ 300 net counts. 
An absorbed power law can be fit also to these data, result¬ 
ing in an absorbing column Ng m (1.7±0.9) x 10 21 cm -2 and 
photon index T = 1.9 ±0.4 (see Fig. [4]). This mean spectrum 
was used to convert the observed count rates of the faintest 
sources into fluxes, see Table[3] 

The number of sources with £^- 5_10keV ) ^ 10 39 erg s -1 
is 8. |Liu| ( |2011 ) classify 9 sources as ULXs, however the two 
lists have only the 5 brightest objects in common (SI, S4, S5, 
S6, Sll). The faintest ones are not in the Liu (2011) list since 
they use a cutoff luminosity of 2. x 10 39 erg s -1 in the 0.3- 
8 keV band and a different distance of D — 36.8 Mpc. The 
other sources in the |Liu (20fT| list which are not listed here 
are: NGC2276-X4, which is the F8 star SAO 1148 (e.g. from 
the Simbad Astronomical Databas^J; NGC2276-X6 which 
is outside the optical region of the galaxy (D 25 = 177V1 from 
NED) at rsj 2' from the nucleus of the galaxy; NGC2276- 
X7 which is on the border of the CCD and did not sur¬ 
vive our thresholds; NGC2276-X9 which is on the next chip 
S2 (CCD6) which we did not even analyze. A comparison 
of XMM-Newton and Chandra images of the field of the 


XMM-Newton ULX candidate XMMU J072649.2±854555 is 
shown in Fig. [5] Instead of a single bright source, Chandra 
found 5 distinct sources (SI, S3, S4, S5, and S6; all in the 
ULX luminosity range, plus a possible other source near SI, 
which is visible also in Fig. [l] but the detection algorithm is 


http: / / simbad.u-strasbg.fr/simbad/ 
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Figure 3. Spectroscopy of the 4 brightest ULXs (SI, S4, S5, S6). 
For each source (as indicated in each plot), we show the spectrum, 
the best-fitting power-law model (solid line) and, in the bottom 
panel, the ratio to the model. 



Energy (keV) 


Figure 4. Combined spectrum of the fainter sources (see Sec¬ 
tion]^]). The solid line indicates the best- fitting power-law model 
and in the bottom panel we plotted the ratio to the model. 



Figure 5. Field of the ULX candidate XMMU J072649.2+854555 
as imaged by Chandra (left) and XMM-Newton (right) in the 
0.3-8 keV energy band. Each image side is 55 X 46 arcsec wide. 
The green contour levels are logarithmically spaced from 0.1 
cts/sq.pix. Several distinct X-ray sources are resolved by Chan¬ 
dra. 


unable to pick it up) roughly following the XMM-Newton 
isophotes (see also Wolter et al. 2011). The luminosity of 
these sources ranges from a few l(Fto a few 10 40 erg s -1 
(Tables [ 2 ] and [3| . 

To investigate possible variability of the Chandra 
sources in the region of XMMU J072649.2+854555, we can 
compare the total luminosity from Davis & Mushotzky] 
(20041 (Z+ -10keV = 1.1 x 10 41 +.7 erg s -1 = 5.7 x 

10 4 '+ 2.9 er g s_1 ) to that of the region in the Chandra 
dataset corresponding to the XMM-Newton point spread 
function (PSF), which we chose as a circle of 30 // ra¬ 
dius centered at the XMM-Newton position. This includes 
sources SI, S3, S4, S5 and S6, plus the unresolved sources 
and diffuse component. We extract a spectrum and fit 
it with the Davis & Mushotzky (2004) best fit model 
(diskbb+mekal) , see their table 1 ) leaving the normaliza¬ 
tion free. We obtain an unabsorbed luminosity £^ 5 - 10keV — 
5.1 x 10 4O di 2 . 9 erg s -1 , which is consistent within uncer¬ 
tainties with the expectation. We make the same com¬ 
parison also for the other two sources detected in both 
the Chandra and the XMM-Newton observation, namely 
source XMMU J072718.8+854636 (which corresponds to 
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source S10), and source XMMU J072816.7+854436 (which 
corresponds to S15). If we consider the model of D avis| 
k Mushotzky| ( |2004| for each of the two sources, we find 
that XMMU J072718.8+85463 is apparently fainter by one 
dex (X^ 5 - 10keV — 5 2 x 10 38 erg s -1 ) with respect to 
the expectation of ^ 5 - 10keV — 4.2 x 10 39 di 5 . 7 erg s -1 = 

2.2 x 10 39 d| 2 .9 erg s -1 ). However if we consider a region of 
30" around S10, we detect in Chandra 158 ± 46 counts, 
from the diffuse contribution of the galaxy besides the 

19.2 net counts from S10: therefore it is probable that 
the XMM-Newton detection was similarly contaminated 
by the plasma component. A different situation is found 
for XMMU J072816.7+854436 which is quite isolated: the 
Chandra luminosity for a power law spectrum is about 
a factor of 1.4 higher (Lx 5_10keV — 7.9 x 10 38 erg s -1 ) 
than expected from XMM-Newton : +^ 5_10keV — 1.1 x 
IQ 39 g ?45 7 erg s -1 = 5.7 x 10 38 d 2 2-9 erg s -1 , which however 
we deem not compelling given the uncertainties in spectral 
shape and the different response matrices of the two instru¬ 
ments. In conclusion, the possible small variation in flux be¬ 
tween the observations, given the different spatial resolution 
of the instruments, cannot be confirmed statistically. 

We note that source S 6 is possibly coincident with a 
triple radio source of total Lsghz = 9.51 x 10 36 erg s 4 , 


which Mezcua et al. (2013) suggest could be the first evi¬ 


dence of an extended jet in a non-nuclear source. The pres¬ 
ence of jets in ULXs has been invoked in a few cases (see 
also Cseh et ah|||2014 ) but the evidence is not firm yet. The 
source has a very blue (hard) color in Fig. [ 2 ] which might 
indicate an excess absorption. This however is not enough 
to imply that the source is a background source since high 
levels of absorption are known to be intrinsic to many ULXs. 


3.1 The X-ray luminosity function 

In order to construct the X-ray luminosity function (XLF) 
of NGC 2276, we have to assess two possible sources of con¬ 
tamination: spurious, unwanted sources which add to the 
XLF, and missing sources due to detection uncertainties. 
We also do not include the central source (S7). To account 
for the possible presence of spurious contamination from 
background sources among the Chandra point sources, we 
randomly excluded one of the three faintest sources, follow¬ 
ing the results of Sect ion[3] for a chance coincidence of 1.4 
source at the faintest fluxes. The contamination becomes 
less than 0.2 sources in the entire galaxy, as defined by 
the D25 diameter, at fluxes Lx ^ 10 “ 14 erg cm -2 s -1 . To 
correct for incompleteness, we base our assessment on fig. 


12 from Kim k Fabbiano (2003) in which the detection 


probabilities as a function of background counts and source 
counts are plotted. The background includes both diffuse 
emission from the galaxy and field background and is esti¬ 
mated in the case of NGC 2276 to be at a level of 0.15 counts 
pixel -1 (by averaging total counts in the galaxy region, up 
to D25, after excluding detected sources). We use the sec¬ 
ond panel, relative to a 2' off-axis, since the whole galaxy 
is within this limit. We apply the interpolated correction to 
all sources with less than 22 counts. This is a conservative 
estimate. In Fig.[ 6 ] we compare the resulting XLF with that 
of the Cartwheel galaxy, one of the richest galaxies in ULXs 
(Wolter k Trinchieri 2004). The Cartwheel XLF was cor¬ 
rected with the same procedure: in this case the background 


-VI-1-1-1- 

v 1 


1.5 - 


\ 

- \ 


j 

zf 

l_J 

o 

3 


0.5 - 


30 


d \ 

\ * ° \ 

\ Os 

\ £ o \ 

\ A O \ 

\ ^ o \ 

\ A A \ 0 O 

\ A \ O 

\ A SO 

\ A \o 
\ A \o 
\ \ 

\ A 'O 

\ ^ 

\ \ 

39 40 

LOG[L(2- 10 keV)] 


Cartwheel 


NGC 2276 


41 


Figure 6. The X-ray luminosity function of NGC 2276 computed 
in the 2-10 keV band by using the sources detected in the galaxy 
area by Chandra. The nucleus has been removed, as well as one of 
the faint sources, to account for the expected contamination from 
background sources (see Sect ion[3|) . The solid lines represent the 
universal XLF for HMXB by |Grimm et al.| |2003), normalized to 
5 (long-dashed line) and 15 Mq yr -i (short-dashed line). 


is lower (0.05 counts pixel -1 ) and therefore we corrected 
only sources with less than 10 counts. Just to give an indica¬ 
tion, the slope of the resulting XLF is c^n = —0.71 + 0.03 for 
NGC 2276 and ac = —0.75 ± 0.05 for the Cartwheel. These 
numbers become lower (qn = — 0.66 and ac = —0.60) if 
we restrict the fit to the log Lx + 39.8, consistent with the 
‘universal’ luminosity function slope of a = —0.60 found 
by e.g. | Grimm et al.||2003| |Swartz et al.|| 201 l] |Mineo et al.| 
2012a for high mass X-ray binaries (HMXBs) and proposed 
to scale with the SFR. The XLF from |Grimm et ah (2003) is 
also plotted in Fig. [ 6 ] for two different values of the SFR of 5 
and 15 Mq yr -1 . Also, by comparing the total X-ray lumi¬ 
nosity or the number of ULXs with the formulas of |Crimm| 
et al. (2003), similar results are obtained: SFR = 5 and 
respectively. These values are in the same 


5.5 Mq yr -1 

range of SFR derived in other wavebands (see Sect. [I]). Also 
eq. (20) from Mineo et al. (2012a) gives a predicted num¬ 


ber of sources with Lx > 10 88 erg s -1 = 3.2 x SFR ^ 16, 
consistent with the number of sources found, taking into ac¬ 
count that at least one is spurious and the incompleteness 
of detection at the limit. 


4 THE DIFFUSE EMISSION COMPONENT 

The nuclear source of NGC 2276 was previously detected 


with ROSAT as variable (Davis et al. 

n 40 d| 2.9 erg s ' 1 


at L 


(0.5-2 keV) 
X 


2 x 10 u d 4 5-7 = 10 « 3 2- 9 erg s -i . Notably, the source was 
not detected in the 2004 XMM-Newton observation, which 
set an upper limit of < 1.2 x 10 39 d | 2 9 erg s -1 

on its emission. However, a correct comparison of the flux 
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Log Radius (arcsec) 

Figure 7. Profile in the direction 280-370°, in red vs. the oppo¬ 
site one: 110-200°, in black. The green dashed line indicates the 
extent of the PSF computed via CHART, normalized to the first 
data point. 


would need a detailed modeling of the surrounding emission, 
given the different resolution of the instruments that have 
observed NGC 2276 through the years. Source S7, coincident 
with the optical nucleus of the galaxy, is extended as we show 
in Fig. [ 7 ] where we compare it with the Chandra PSF derived 
with CHART with standard procedures at the position of the 
nucleus. We plot two different profiles in the NE and SW 
directions (110 to 200° and 280 to 370°, respectively) to 
show the different extent at large radii in the two directions, 
as already noted by Rasmussen et al. ( 2006| . 

Since the diffuse emission of NGC 2276 and in its prox¬ 
imity was studied in detail by Rasmussen et al. (2006), we 
just verify the consistence of our results for the diffuse emis¬ 
sion of the disc (after removing the point sources) for their 
region A (see their figure 2). We used an extraction region 
as similar as possible to that of Rasmussen et al. (2006): a 
circle of 50 /x radius, visible in Fig. [l] 

We make a more general hypothesis that there might 
be intrinsic absorption in the galaxy, and we bin the spec¬ 
trum so as to have a minimum of 100 counts per bin, 
but we checked that the values are not very dependent 
on binning. In this way, in addition to the Galactic ab¬ 
sorption, we measured an absorbing column equivalent to 
Ah = ( 4 . 2 ^ 3 ) X 10 21 cm -2 for the best fit (xt — 1.01 
for 13 degrees of freedom (dof)). The resulting spectrum 
is plotted in Fig.[8] This results in a lower temperature of 
kT = 0.18^q'q 2 keV and an only poorly constrained pho¬ 
ton index T = 4 d= 2, but we note that the power-law com¬ 
ponent is necessary to properly fit the spectrum and we 
freeze it to T = 1.7, appropriate to represent both a low lu¬ 
minosity active galactic nucleus (AGN) component and/or 
the unresolved point source component. The corresponding 
total luminosity is £^' 5-2keV ) = 1.8 X 10 41 erg s -1 . The 
luminosity of the power-law component is £^' 5-2keV ) = 



Energy (keV) 


Figure 8. Spectrum of the central region of the galaxy (see Sec¬ 
tion^). The solid line indicates the best-fitting power-law model 
and in the bottom panel we plotted the ratio to the model. 


2.1 x 10 40 ergs -1 , about 12 per cent of the total lumi¬ 
nosity and consistent with the expectation for the unre¬ 
solved point source component derived from the K band 
luminosity of NGC 2276 Lk = 10.7, which corresp onds to 
L 


-2 — 10 keV 

12L- 


10 * 


,40 erg s 1 


by using the relation in 


Kim & 


Fabbiano (2004). Since there should also be a fraction of 


unresolved HMXBs, this leaves little room for a powerful 
AGN. 

We observe that our revised temperature and luminos¬ 
ity have essentially no impact on the results of [R asmussen] 
et al. ( 2006| (actually, the higher luminosity and lower tem¬ 


perature enforce their conclusion that diffuse X-ray emis¬ 
sion from the western edge does not arise from a shock- 
compressed gas component). 

If we compute the total emission from the D25 of the 
galaxy, then the diffuse component is £^- 5_2keV ^ 2.5 x 
-1 : we can compare it with the findings of 


10 


er g s 


|et al. (2012b) that L% ttuse /SFR = 1.5 x 10*~ erg s 


Mineo 
- 1 / 


(Mq yr 1 ), albeit with a large scatter of ~ a dex and very 
model dependent. Again this gives a SFR in the range 5-25 
M© yr -1 . 

We recall also that the source is immersed in an IGM 
component: we report the parameters found by |Davis et al.| 
(1996) using ROSAT data, rescaled to our distance, which 
we will exploit for the subsequent simulations. They found 
an extent of about 25 x (0.24 Mpc) for the X-ray emitting 
IGM, centered close to, but not exactly on, the central 
elliptical galaxy NGC 2300. The surface brightness of the 
medium is described by a King function with core radius 
ro ~ 4.3 / and index (3 = 0.41, and has a luminosity of 
Lx o1 = 1.12 x 10 42 d 2 5 . 7 = 5.8 x 10 41 d 2 2 .9 erg s -1 . The mass 
of the hot intragroup gas is M gas = (1.25 ±0.13) x 10 12 Mq. 


Rasmussen et al. (2006) find similar values, but the ROSAT 


larger field of view is better suited to measure this large 
extent component. 


5 WHERE DOES THE MORPHOLOGY OF 
NGC 2276 COME FROM? 

Whether the perturbed morphology of NGC 2276 is the re¬ 
sult of tidal interactions with NGC 2300, or of ram pres¬ 
sure stripping is still an open question. The simplest way 
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to estimate the importance of tidal forces, without running 
a simulation, is to compare the centripetal acceleration of 
NGC2276 (a c ) to the radial acceleration (a r ) exerted by the 
other galaxies in the group, which is dominated by NGC2300 
and by the hot gas component: 


a r (D) — G Mgroup(D) 


a c (r ) — 

1 

D 2 


GM(r) 


1 


(D + r) 2 


(1) 

( 2 ) 


where r and M(r) are the three-dimensional distance from 
the center of NGC2276 and the mass of NGC2276 inside 
r, respectively, while M group (D) is the mass of the group 
enclosed within a distance D from the center of NGC 2300. 

To derive the actual distance between NGC 2276 and 
NGC 2300 we used the velocity of NGC 2276 with respect 
to the IGM, derived assuming that the IGM is coupled to 


NGC 2300. Rasmussen et al. (2006) derive a Mach number 


of M — 1.70 ± 0.23, corresponding to a three-dimensional 
velocity of v = v s X M — 865 db 120 km s -1 . The distance 
D = x/cosa is then derived by measuring x = 83.2 kpc (the 
distance between the centers of the two galaxies projected 
on the plane of the sky) and cosa = Av r /v , where Av r = 
421 km s _1 is the relative recession velocity between the 


two galaxies (Ras mussen et al.|2 006). The three-dimensional 
distance between the two galaxies is then 170 kpc, under 
these assumptions. 

For r — 13 kpc (which is approximately the isopho- 
tal radius of NGC 2276) and M(13kpc) = 2 x 10 11 Mq, 
we obtain a c — 1.6 X 10 -8 cm s -2 . For D = 170 kpc, 
a r (D) ~ 5 x 10 10 cm s 2 . Thus, we expect tidal forces to 
be small today, when compared with the self-gravity of the 
galaxy. Let us consider now the importance of ram-pressure 
stripping. The distance from the center of a galaxy at which 
ram-pressure stripping becomes important can be evaluated 
as ( Yoshida et al.| 2008) 


_ h ^ / G M* Mism 

2 U \2 7 TV 2 plGM h 4 


(3) 


where h is the scale-length of the disc, M* and Mism 
are the total stellar mass and inter-stellar medium mass of 
the galaxy, v is the relative velocity between the galaxy and 
the IGM, and pigm is the local density of IGM. If we as¬ 
sume h — 3 kpc, M* — 2.8 x 10 10 Mq, Mism = 8 x 10 9 Mq, 
v = 850 km s _1 and pigm = 5 X 10 -27 g cm -3 , which are 


consistent with the observations of NGC 2276 (Davis et al. 
Ras mussen et ak||2006 ) and its environment (we use 


1996 


the same values for the simulations described in the next 
Section), we obtain r rp = 6.5 kpc, which is larger than h 
but still inside the isophotal radius of NGC 2276 (13.2 kpc). 
Thus, we expect ram-pressure to be quite effective in per¬ 
turbing the morphology of NGC 2276. 

Even if from equations 0 and 0 we have shown that 
tidal forces are not as efficient as ram pressure, it is reason¬ 
able to expect that there is an interplay between tidal forces 
and ram pressure, leading to an enhancement of NGC 2276 
stripping with respect to the case in which ram pressure was 
the only process at work (see e.g. the argument in Yoshida 
et al.|2008 ). In addition, viscous stripping (e.g. Nulsen|1 982) 


might also play a role, removing matter from the ISM. 


6 NUMERICAL SIMULATIONS 

To study the morphology of NGC 2276 in more detail 
than with our order-of-magnitude calculations, we ran N- 
body/smoothed-particle hydrodynamics (SPH) numerical 
simulations of the interaction between NGC 2276 and its 
environment. We used the gadget-2 code (Springel et al. 
2001 Springel||2005 ). 


The initial conditions are taken to the simplest possible 
level. In all the simulations, we put NGC 2276 on a parabolic 
orbit around NGC 2300 (assumed to sit at the bottom of the 
gravitational potential well of the group). The two galaxies 
start from an initial separation of 170 kpc, 85 Myr before 
the periapsis. At periapsis, the galaxies reach their minimum 
distance of 160 kpc, and after another 85 Myr they reach 
their present separation of 170 kpc. We ran the simulations 
for 1 Gyr, but the most relevant effects occur during (and 
within < 200 Myr after) the periapsis passage. 

We describe the elliptical galaxy NGC 2300 as a rigid 
potential, to reduce the computing time. We further assume 
that the total contribution is due to the dominant compo¬ 
nent of dark matter, described by a Navarro, Frenk and 
White (NFW) profile ( Navarro et al.| 1997): 

Po 


p(x) = 


(4) 


X (1 + x) 2 ’ 

where x = M with r s = R 200 /C and R200 = (M to t G/Hq) 1/3 
(Ho — 73 km s -1 Mpc -1 being the Hubble constant, M to t 
being the total mass and c being the concentration param¬ 
eter), and 

_ 200/9 c ritC 3 , . 

P0 ~ 3 (ln(l + c)-c/(l + c))’ l j 

where p C rit is the current critical density of the Universe. We 
assume M to t = 1.3 x 10 13 Mq and concentration parameter 
c = 10. 

The second ingredient is the IGM. We perform three 
different simulations for different choices of the IGM profiles. 
The first is a ‘dry’ simulation, with no IGM. This is a control 
simulation aiming to isolate the purely tidal effects between 
the galaxies. 

The second and third simulations assume that the IGM 
is present, but with different distributions. In both cases 
the distributions are centered on NGC 2300. In both simula¬ 
tions, the gas is assumed to be adiabatic, with initial temper¬ 
ature 2 x 10 6 K, the total mass of gas is M g = 1.25 x 10 12 Mq 
and the mass of a single gas particle is 10 6 Mq (for a total 
number of 1.25 x 10 6 particles). 

In the second simulation, the IGM has a NFW profikQ 
with concentration c = 10. 

Finally, in the third simulation the gas is distributed as 
a King model ( King|p~ 962): 

-3 /3+0.5 


p(r) = po 


1 + - 


(A 


( 6 ) 


7 The NFW profile usually refers to dark matter, and not to 
baryonic matter. The rationale of adopting a NFW profile for the 
baryonic IGM is the following. For large radii, the NFW has the 
asymptotic form g oc r ~ 3 to be compared with the g oc r~ 2 of 
an isothermal sphere. The NFW would model an IGM distribu¬ 
tion concentrated close to NGC 2300, while the isothermal sphere 
models a more widespread IGM. 
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where po — 3 .1 x 1 0 6 cm 6 (from Davis et al. 1996), 


/3 = 0.41 (from Rasmussen et al.||2006 ) and r c = 4.0' (from 
Ras mussen et al. 2006). 

The third ingredient is a model of NGC 2276. We 
parametrize it, according to Widrow & Dubinski (2005), as 


the sum of three different components: bulge, disc and dark 
matter halo. The total number of particles in the dark mat¬ 
ter halo, bulge and disc of the NGC 2276 model are 2 x 10 6 , 
10 5 and 6 x 10 5 , respectively. The softening length is 100 pc. 

The dark matter halo is modeled as a NFW potential 
with M to t = 6 x 10 11 Mq and c = 10. The bulge is rep¬ 
resented by_a Sersic profile, according to the Prugn iel fe| 


Simien 


(1997) deprojection pb(r) oc (r/r e ) p exp (r/r e ) 1/n , 


where we assume r e — 0.83 pc, p — —1 and n — 1.83. The 
normalization of the profile is set by the total mass of the 
bulge Mb = 6 x 10 9 Mq. 

The gas and the stellar component are distributed as 
an exponential Hernquist disc ( |Hernquist|1 993): 


Pd(R,z) = 


M d 

47vh 2 zo 


exp ( s—R/h ) sech 2 (z/zq) . 


(7) 


where R and z are the cylindrical coordinates, M d = 3.6 x 
10 10 Mq is the total mass of the disc, h — 3 kpc is the scale- 
length of the disc and zo = 0.3 pc is the scale-height of the 
disc. The disc is composed of both gas and stellar particles, 
with the total mass of gas and stars being Md, g = 0.8 x 10 10 
and Md,* = 2.2 x 10 10 Mq, respectively. 

The gas component of the NGC 2276 model is assumed 
to be adiabatic, with initial temperature of T = 2 x 10 4 K 


(Rasmussen et al. 2006). We do not introduce a cooling 


function, but verify ‘a posteriori’ if the cooling time t C ooi 
is shorter than the dynamical free fall time tdyn to estimate 
the rate of star formation. In particular, we estimate the 
cooling time as 


Cool - p. 


nkT 


2 n e n p A (Z, T, n) ’ 


( 8 ) 


where n is the total gas density, n e is the density of electrons, 
n p is the density of protons, k is the Boltzmann constant, 
T is the temperature, and A(Z, T, n) is the cooling function, 
as a function of n, T and of the metallicity Z. We used the 
cooling function of Sutherland & Dopita (1993). 


We assume that gas is able to form stars if 


Cool ^ ffiyn ^ t\ ife? (9) 

where tn f e is the life time of a gas clump or structure, derived 
directly from the simulations. 

Here below we shortly summarize the results of our sim¬ 
ulations. 


No gas is stripped from the galactic disc, and the tails 
observed by Rasmussen et al. (2006) are not reproduced. 


6.2 NFW simulation 

The evolution of the disc galaxy changes significantly, if 
the IGM is accounted for. As we expect from our order- 
of-magnitude estimate, ram-pressure stripping and viscous 
stripping are efficient, and lead to the formation of tails al¬ 
ready during the periapsis passage (see Fig.|9|. We also see 
that the front side of the galaxy disc is perturbed, recalling 
the asymmetric ‘bow-shock-like’ shape of NGC 2276. 

The total mass in the tails is ~6 x 10 8 Mq after the 
first periapsis passage (it will rise to ~1 x 10 9 Mq at the 
end of the simulation). A part of this mass forms clumps 
that are visible in Fig. [9] The cooling time in the galaxy 
disc is ~10 4 - 5 yr, while it is 10-20 Myr in the ram-pressure 
stripped tails. The dynamical time tdyn is 200 Myr and 10 
Myr in the galaxy disc and in the tails, respectively. Thus, 
star formation occurs in the galaxy disc and might occur 
also in the tails. 

Finally, a small but significant fraction of the gas com¬ 
ponent of the NGC 2276 model (8.6 x 10 s Mq, i.e. 6.6 per 
cent of the entire gas component) is completely stripped and 
lost from the galaxy during the entire simulation. 


6.3 Isothermal sphere simulation 

The third simulation is very similar to the second one. Even 
in this simulation, ram pressure produces a deformation in 
the front side of the galaxy disc and leads to the formation 
of tails and clumps. The total mass in the tails is 6 x 10 8 
Mq after the periapsis passage, and raises to 7.5 x 10 8 Mq 
at the end of the simulation, with no significant difference 
with respect to the second simulation. In the tails, t coo l ~ 
tdyn ~ 10 Myr, suggesting that star formation is possible. 
In the disc galaxy, Cool ~ 10 4-5 yr << tdyn ~ 300 Myr, 
indicating that star formation should occur. 

The total mass of gas that is lost from the disc galaxy by 
the end of the simulation is 1.5 x 10 8 Mq, i.e. ~ 12 per cent 
of the entire initial gas mass in the NGC 2276 model. This is 
nearly a factor of two larger than in the second simulation. 
The difference can be explained with the fact that the IGM 
density goes as pigm oc r~ 3 and oc r ~ 2 in the second and in 
the third simulation, respectively. Since the density of the 
IGM at large distances from the elliptical galaxy in the third 
simulation is higher than that in the second simulation, we 
expect that ram pressure stripping will last longer than in 
the second simulation (i.e. it is important even when the 
disc galaxy is far from periapsis). 


6.1 Dry simulation 

Our order-of-magnitude estimate shows that tidal forces 
should not be very important for NGC 2276. On the other 
hand, even in the simulation without IGM we observe the 
formation of (not strongly pronounced) tidal arms in both 
the stellar and the gaseous components. In addition, the gas 
disc thickens. The regions where condition [9] is satisfied are 
very rare in the dry simulation, and we expect that star for¬ 
mation is not enhanced by the interaction, even during the 
periapsis passage. 


6.4 Summary of simulation results and caveats 

In conclusion, we found that ram-pressure stripping and 
viscous stripping are the dominant processes and can lead 
to both the formation of the tails and the deformation of 
the front side of the disc galaxy. Tidal forces contribute 
marginally, by producing tidal arms and thickening the 
gaseous disc. We found a factor of 2-5 lower SFR than in the 
observations, but we are aware that our very approximate 
treatment of gas underestimates both cooling and SFR. 
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Figure 9. Color-coded density map stars (top panel) and gas 
(central and bottom panel) at t = 300 Myr in the NFW simulation. 
The boxes measure 300 kpc per edge. Top panel: the density of 
stars has been projected along the z-axis; central ( bottom) panel: 
the density of gas has been projected along the z-axis (y-axis). A 
tail structure is well visible in the central and bottom panel. 


Moreover, the version of the SPH technique adopted in 
GADGET- 2 is known to poorly resolve dynamical instabilities 
and mixing (e.g. Agertz et al. 2007). The main reason is 
that the SPH formulation adopted in GADGET- 2 inaccurately 
handles the case of low-density regions in contact with high- 
density ones. As a consequence, gas turbulence is artificially 
suppressed, and gas stripping due to ram pressure might be 
underestimated with respect to other formalisms (e.g. grid 
codes). In particular, we expect that the gas mass in the 
stripped tails is underestimated in our simulations, and the 
morphology of the resulting tails is quite smoother than it 
should be. This implies that our results can be regarded 
as lower limits for the ram-pressure effects. Bearing these 
caveats in mind, we conclude that our simulations are able 
to qualitatively reproduce the peculiar shape of NGC 2276, 
and to confirm the importance of ram pressure and viscous 
stripping. 


7 SUMMARY & CONCLUSIONS 


We have analyzed a Chandra observation of NGC 2276, a 
spiral galaxy in the group of the elliptical NGC 2300. We 
have paid particular attention to the point sources detected 
in the area of the galaxy. We find that a large number of 
ULXs are present, which makes NGC 2276 one of the richest 
environments for this still enigmatic bright sources. 

We have seen in W olter et al.| (2011) that six of the X- 
ray bright sources, four of which are ULX s, are positionally 
coincident (l // -2 // ) with Hll regions (from Hodge & Kenni- 
cutt|1983]|Davis et al.| 1997). However the spatial resolution 


at the distance of NGC 2276 is not enough to warrant a 
physical association. At least five supernovae have been ob¬ 
served in the last fifty years (see Section[l] for references), 
two close to the center and two in the western region where 
ULXs are found. 

All observations concur on linking the high activity 
with the presence of a large number of ULXs, although 
there is no one-to-one correspondence between peaks in dif¬ 
ferent bands, as is already seen in other galaxies, like the 
Cartwheel ( Wolter fe Trinchieri|2004 ). Radio emission is en¬ 
hanced throughout the galaxy body: one of the sULX is 
possibly associated to either a radio bubble or an extended 
jet structure ( Mezcua et al.|2 013), however the complex ra¬ 
dio field makes the interpretation difficult. The source has a 
very hard spectrum: we cannot exclude the possibility that 
it is a chance superposition with a background AGN. 

We detect 15 sources associated with the galaxy 
NGC 2276 in the Chandra observation, out of which 8 are 
above the ULX luminosity threshold. We derive an XLF 
and confirm the previous estimates of the SFR of NGC 2276 
in the range SFR = 5-15 Mq yr -1 . The value is close to 
those of the Antennae and t he Cartwheel, which are fitted 
with SFR = 7.1 yr -1 flZezas &; Fabbia no 2002) and 
~ 20 Mq yr -1 (Wolter & Trinchieri 2004), respectively. 
These numbers reconcile expectation from Mapelli et al.| 


(2010) with observations. 


We know that moderate variability is rather common 
in ULXs (see e.g. Crivellari, Wolter & Trinchieri 2009] and 
Zezas et al. 2006) and in recent years a few transient ULXs 
have also been identified (e.g. Sivakoff et al.||2008 Middle- 


ton et al.|2012 2013 Soria et al.|2012 ; Esposi to et al.| 20*13 ). 
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We investigate, as much as possible given the different in¬ 
struments used with widely different spatial resolution and 
response, flux variations for the bright source detected in 
XMM-Newton. We do not find any significant evidence for 
variablity, although a better comparison could be done by 
using the same instrument repeatedly. 

We confirm parameters for the diffuse emission as de¬ 
rived by Rasmussen et al. (2006), however we find that 
the spectrum is best fitted by adding an intrinsic absorp¬ 
tion component that raises the unabsorbed luminosity to 
^( 0 . 5-2 keV) _ ^ g x iq4i er g g -i p ea k of this distribu¬ 

tion is at the center of the galaxy. The spectrum is consistent 
with the sum of a hot plasma diffuse component, plus an un¬ 
resolved component due to active stars and binaries account 
for a small fraction of the total luminosity (see Section^. 

The Chandra data display no point source at the cen¬ 
ter of the galaxy. A low luminosity variable AGN of Lx ^ 
few x 10 39 erg s -1 might lurk in the midst. However the ob¬ 
servational results of both the ROSAT HRI (Total Lx == 


Davis et al. 


1997) 


3.2x10 tt 4 5 . 7 erg s in a 22 radius circle: 
and the XMM-Newton (Lx ^ 2.4 x 10 39 d 4 5-7 erg s _i above 
the diffuse emission: Davis X Mushotzky|2004 ) datasets are 
consistent with no point source in the center of the galaxy. If 
anything, the presence of two ULX (Sll and S12) at r ^ 30 // 
could be responsible for some variability on long time scales. 

We have performed a number of simulations to derive 
the source of activity and asymmetry in NGC 2276 and con¬ 
firm that this is linked to the large number of ULXs found. 
Assuming a simplified scenario we derived a morphology 
very reminiscent of the observed one and we derive a total 
mass within a factor of 2 from observations. The presence of 
a gas envelope around the central elliptical galaxy NGC 2300 
is necessary in order to produce the effects, along the same 
lines as suggested by Rasmussen et al. (2006). 

In conclusion, we found that ram-pressure stripping 
and viscous stripping are the dominant processes and can 
lead both to the formation of the tails and to the deforma¬ 
tion of the front side of the disc galaxy. Tidal forces con¬ 
tribute marginally, by producing tidal arms and thickening 
the gaseous disc. Therefore, even if at a simple level, we con¬ 
clude that our simulations are able to reproduce the very 
peculiar shape of NGC 2276 and the mass of the tidal tails. 
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